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Abstract 



We study the masses and decay constants of heavy-light flavour mesons D, D s , and B s 
in a QCD Potential model . We use the mesonic wavefunction in momentum space and estimate 
the pseudoscalar decay constants. We also calculate leptonic decay widths of these mesons to 
compute the leptonic branching ratio for different leptonic channels. The results are found to be 
compatible with available data. 
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1 Introduction 

The study of the wave functions of heavy-flavored mesons like B and D are important not only for 
studying the properties of strong interaction between heavy and light quarks, but also for investigating 
the mechanism of heavy meson decays. The wave function determines the momentum distributions 
\Tl \ of the quark and antiquark in mesons, which is an important quantity for calculating the amplitude 
5_i ; of heavy meson decays [1, 2]. Heavy hadron spectroscopy has played a major role in the foundation 
— ■ of QCD. In the last few years however it has sparked a renewal of interest due to the numerous 
data available from the B factories, CLEO, the Tevatron and by the progress made in the theoretical 
methods. The remarkable progress at the experimental side, with various high energy machines such 
as BaBar, BELLE, B-factories, Tevatron, ARGUS collaborations, CLEO, CDF etc., for the study of 
hadrons has opened up new challenges in the theoretical understanding of light-heavy flavour hadrons. 



Regarding the value of the Pseudoscalar decay constants, experiments and 1QCD calculations agree 
very well with each other on the value of /o.But the most precise Nf = 2 + 1 lattice QCD (1QCD) 
result from the HPQCD/UKQCD collaboration [4]for f Ds = 241 ± 3 MeV was found to be about 3a 
less than the PDG average(2008) value for fo s which is 273±10Me\^[3] . The discrepancy concerning 
fn s is quite puzzling because whatever systematic errors have affected the 1QCD calculation of fjj , 
should also be expected for the calculation of fc s [5]. However, the discrepancy is reduced to 2.4a 
with the new (updated) data from CLEO [6,7] and BaBar [8], together with the Belle measurement 
[9] and the latest PDG average is fjj a = 257.5 ±6.1 MeV [10]. Lately, the HPQCD collaboration has 
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also reported its updated result as fo s = 248.0 ± 2.5 MeV[ll]. 



Though there exist many potential models with relativistic and non relativistic considerations em- 
ployed to study the hadron properties based on its quark structure [12-14], the most commonly used 
potential is the combination of coulomb term and linear confining term V(r) = :z ^ s - + br + c [15,16,17]. 
The present authors have been persuing a specific potential model considering this potential with its 
Coulombic part as perturbation[18] as well as linear part as perturbation [17,19,20]. In this work, 
we have transformed the wavefunction with linear part as perturbation in momomentum space by 
applying Fourier transformation. This wavefunction is then used to study the decay constants and 
leptonic branching ratio of B and D mesons in this QCD potential model approach. 

We discuss the formalism in section 2 and summarise the results and conclusion in section 3. 



2 Formalism 

2.1 Wave function in the model 

We start with the ground state (/ = 0) spin independent non relativistic Fermi-Breit Hamiltonian 
(without the contact term) 

H ---,~ + br + c - (1) 

With the linear term br + c as perturbation and using Dalgarno method,the wave function in the 
model is obtained as [19,20] : 

N> = 2 * (3) 

(3 - 26) C* - Ip&agr (5 - 2c) C + ^W T (7 - 2e)) 



C' = l + cA ^ixal (4) 

^ = ~ I 1 ( 5 ) 



rrii + m 
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a = ( -flag ) (6) 



4 ^ 2 



The QCD potential is taken as 



e=l-\ll-[-a s ) . (7) 



V(r) = -^a s + br + c (8) 
3r 
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Here A Q is the undetermined factor appearing in the series solution of the Schrddinger equation. However 
with A =0 the effect of c vanishes entirely from the solution. 

The wavefunction in momentum space can be obtained by using the Fourier transform as 

i> (P) = 1 3/2 / d z re=%ri) (r) . (9) 

Separating the variable-dependence of the momentum space wave function as 

^(p)=^(p)Y lm (9,<P) (10) 
one can obtain for / = in the natural unit as[16] 



i> (p) = y J drsin (P r ) ^( r )- (11) 
Then using Eqs.(2)and(ll) and the standard result 

J^-m™**-^^, (12) 



one can obtain the normalised wavefunction in momentum space as 

c 



( . iVv / 2(2-e)r(2-e) 

Wrel+conf [P) = 



(4 - e) (3 - e) fiba 3 



2 (1 + a 2 p 2 ^ 



(13) 



This simplified form of the wavefunction gives the momentum distribution of the quark and anti 
quark. 

2.2 Masses and Decay constants of D and B mesons 

In the relativistic quark model, the decay constant can be expressed through the meson wave function 
ipp(p) in the momentum space [21,22]as 



rnr r_fi 

V MpJ (2tt 



with Xp = — 1 for pseudoscalar mesons and E q = \/p 2 + q 2 . 

The pseudoscalar mass M P of the mesons are calculated by using the relation [23, 24] 

M p = itlq + m-Q + (H) (15) 
where the expectation value of the hamiltonian is 

(H) = (£-) + (V(r)). (16) 

The strong running coupling constant appeared in the potential V(r) in turn is related to the 
quark mass parameter as [19, 25] 

«> M = 7 ^ , MU ,s ( 17 > 



(ll-J^ta^) 



mesons 


present 


experimental masses jlOj 


D(cu/cd) 


1870.82 


1869.6 ±0.16 


D(cs) 


1966.62 


1968 ±0.33 


B d (bd) 


5273.50 


5279 ± 0.29 


B s {bs) 


5365.99 


5366 ±0.6 



Table 1: Masses of heavy-light mesons in this work and comparison with experimental data. All 
values are in units of MeV. 





fD 


Id s 


Sb 




Present 


205.14 


241.84 


201.09 


292.04 


Experiment [28,29] 


206 ±8.9 


260.0 ±5.4 


204 ± 31 




LQCD [30] 


218.9 ± 11.3 


260.1 ± 10.8 


196.9 ±8.9 


242 ± 9.5 


LQCD [11] 


213 ±4 


248 ± 2.5 






ExChQm[27] 


207.53 


262.56 


208.13 


262.39 


LC [28] 


206 ±8.9 


267.4 ± 17.9 


204 ± 31 


281 ± 54 


LQM [31] 


211 


248 


189 


234 


FC [32] 


210 ±10 


260 ± 10 


182 ±8 


216 ±8 


BS [33,34] 


230 ± 25 


248 ± 27 


196 ± 29 


216 ±32 


RQM [22] 


234 


268 


189 


218 


RPM [35] 


208 ± 21 


256 ± 26 


198 ± 14 


237 ± 17 



Table 2: Decay constants of pseudoscalar heavy-light mesons(in MeV) computed in this work 
and comparison with experimental [28, 29] and theoretical results from (2+l)-flavour asqdat ac- 
tion[30],HPQCD[ll], extended chiral quark model(ExChQm) [27], Light cone wavefunction[28], light- 
front quark model (LQM) [31], field-correlator method (FC) [32], Bethe-Salpeter method (BS) [331134] , 
relativistic quark model (RQM) [22] ,relativistic potential model(RPM)[35] 

where, n/ is the number of flavours, fi is the renormalisation scale related to the constituent quark 
masses as a = 2 m ', mj . Mr is the background mass related to the confinement term of the 

potential as Mb = 2.24 x 5 1 / 2 = 0.95GeV. In this calculation we have taken nf = 3 as in 
ref. [25,26] and Aqcd = 0.200GeV . The input parameters used in the numerical calculation are 
m d = 0.336GeV, m s = 0.465GW, m c = 1.55GeV, m b = 4.97GeV with b = 0.183GW 2 , c=-0.37 GeV 
and cA =lGeV 2/3 . 

Using eq.15 and eq.16 we compute the pseudoscalar ground state masses of the heavy light pseu- 
doscalar mesons and compare with the experimental data in Table. 1. The results are found to be in 
good agreement with the experimental data. Again, using these computed masses we employ eq.14 
to obtain the pseudoscalar decay constants. The results are then compared with the available ex- 
perimental and theoretical values in Table. 2. The results are found to be compatible with available 
experimental and theoretical values. 

We note that the present result, ExchQm [27] and that from LC [28] give fs s > 260 MeV, while 
other results give fs s < 240 MeV. Hence, the experimental measurements for fs s can be a good 
testing ground for theoretical reliability of each model as shown here. 
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mesons 


BRr x 10~ 3 


BR U x 10~ 4 


BR e x 10" 6 


D(cd) 
Expt.[10] 


1.08 [present] 
< 1.2 


3.89 [present] 
3.82 ±0.32 ±0.09 


0.092 [present] 
< 8.8 


B. Patel etal.,[36] 


0.9 


6.6 


0.015 




BRr x 10~ 2 


BR U x 10~ 3 


BR e x 10~ 4 


D(cs) 
HFAG [37] 


5.43 [present] 
5.38 ±0.32 


5.33 [present] 
5.8 ±0.43 


0.00 13 [present] 


Expt.[10] 


5.6 ± 0.4 


5.8 ± 0.4 


< 1.2 


B.Patel etal.,[36] 


8.4 


7.7 


0.0018 




BRr x 10~ 4 


BR^ x 10" 6 


BR e x 10~ 6 


B(bu) 
Wolfgang etal.,[38] 


1.07 [present] 
0.80 ±0.12 


0.48 [present] 


0.0001 [present] 


Expt.[10] 


1.8 ±0.5 


< 1.0 


< 1.9 



Table 3: Leptonic branching ratio of D, D s and Bd mesons for three leptonic channels and comparison 
with experiment and theoretical results. 

2.3 Leptonic decay width and Branching ratio of D, D s and B mesons 

Charged mesons formed from a quark and anti-quark can decay to a charged lepton pair when these 
objects annihilate via a virtual W ± boson. Quark-antiquark annihilations via a virtual to 
the final states occur for the D^ 1 and B^ mesons. Purely leptonic decay processes are rare 

but they have clear experimental signatures due to the presence of a highly energetic lepton in the 
final state. The theoretical predictions are very clean due to the absence of hadrons in the final state. 
The partial decay width for P — >■ iv reads: 

Y{P^iv) = ^flM pm \ (l-^)V//, (18) 

where Gp, P, fp, M P , and denote the Fermi constant, generic pseudoscalar(PS)meson, PS- 
meson weak-decay constant, PS-meson mass and lepton mass respectively. Vf g stands for the CKM 
matrix element for the quark flavors / and g. The importance of measuring r (P — > Iv) depends 
on the particle being considered. In the case of the B~ meson, the measurement of V (B~ — > t~v) 
provides an indirect determination of |T4&| provided fp is given by theory. 

The leptonic widths of the charged PS mesons are obtained by using eq.18 and employing the 
predicted values of the pseudoscalar masses and decay constants fp, fp g and fp . The leptonic 
widths for separate lepton channels by the choice of m; =TjAt e are computed to obtain the branching 
ratio of D and B mesons. The branching ratio of the heavy-light mesons are calculated by using the 
relation 

B = T P T(P ->■ tv). (19) 

The life time of these mesons are t d = 1.04ps, t Ds = 0.5ps, r B = 1.63ps and the CKM elements 
Vcd = 0.230, V cs = 1.023, V u b = 3.89 x 10~ 3 are taken from the world average value reported by 
Particle data group [10]. The present results as tabulated in Table. 3 are in accordance with the 
available experimental values. 
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3 Summary and Conclusion 



In this work, we have computed the Pseudoscalar masses and decay constants of heavy-light mesons(B 
and D). We have transformed the wavefunction from r space into momentum space and have used 
it to obtain the weak decay constants with its relativistic effect. These masses and decay constants 
are then used to compute the branching ratio of D, D s and B mesons for the three leptonic channels 
r, /j, and e. B s meson being neutral in nature, does not show leptonic decay and hence the leptonic 
branching ratio for B s meson is not computed. 



• The renormalization scale of the model was set to be Kqcd — 200 MeV, with the approximation 
that the Aqcd for the heavy-quark effective mass is the same as that for the light quark, taking 
into account that the QCD-vacuum structure is not affected much by the heavy quarks, although 
the heavy sources may distort the vacuum to a certain extent. 

• The ground state masses of D and B mesons computed in this approach are found to be well 
consistent with the experimental vaues. 

• We obtain the decay constants as f D ,D s ,B,B s = (205.14, 241.84, 201.09, 292.04) MeV which are 
qualitatively compatible with available experimental and theoretical values. Except Jb s , other 
values of the decay constants locate inside the experimental errors. However, with a variation 
of A.qcd for D and B mesons one obtains more compatible results with the data, although we 
do not provide those numerical numbers. 

• The computed value in the present work = 1.178 is found to be in good agreement within 

the error limit of the recent result Lattice(HPQCD) f -f^ = 1.164 ± 0.018[11] and Lattice(FNAL 

and MILC)^ = 1.188 ± 0.025[39]. However the result of ^ = 1.45 are found to be larger 
than the other theoretical values. 

• The leptonic branching ratio calculated in the present work for three leptonic channels are 
comparable with their empirical and PDG average data. The large experimental uncertainity 
in the electron channel make it difficult for any reasonable conclusion. Furthermore, the ratio 
of branching ratio in the present work is found to be R = f !^1^ T , u \ = 10.18 which is not far 
from the experimental result 9.2 ± 0.7 and Standard Model result 9. 76 [29]. 

Taking into account all the results summarized above, we can conclude that the present theoretical 
framework of Potential Model is a qualitatively successful model to study the heavy-light Pseudoscalar 
mesons. From a phenomenological point of view, the present theoretical framework is a considerably 
useful tool to investigate various physical quantities for the heavy-light quark systems, such as the 
Isgur-Wise function, heavy-light meson coupling const ants, form factors and Charge radii and so on. 
Related works are under progress and will appear elsewhere. 
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